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We evaluate new observational constraints on the two-parameter scale-dependent spectral index 
predicted by the running-mass inflation model by combining the latest Cosmic Microwave Back- 
ground (CMB) anisotropy measurements with the recent 2dFGRS data on the matter power spec- 
trum, with Lyman a forest data and finally with theoretical constraints on the reionization redshift. 
We find that present data still allow significant scale-dependence of n, which occurs in a physically 
reasonable regime of parameter space. 



I. INTRODUCTION. 

The recent observation of the cosmic microwave 
background (CMB) angular anisotropics [|, |[ ||, f|, 
0, 0, 0| have revealed an outstanding agreement be- 
tween the data and the inflationary predictions of 
a flat universe and of a primordial scale invariant 
spectrum of adiabatic density perturbations (see e.g. 

Furthermore, CMB measurements are in perfect 
agreement with other types of data and a global con- 
sistent picture is slowly emerging. For example, the 
CMB constraint on the amount of matter density in 
baryons, u>b, is now consistent with the independent 
constraints from standard big bang nucleosynthesis 
(BBN) (see e.g. Q). Also, the detection of power 
around the expected third peak, on arc-minutes scales, 
provides a new and independent evidence for the pres- 
ence of non-baryonic dark matter (see e.g. jlaj). At 
the same time, early data releases from the 2dFGRS 
and SDSS galaxy redshift surveys Jl^, |l9| are liv- 
ing up to expectations and combined analysis of all 
these datasets are placing strong constraints on most 
of the cosmological parameters ^0|, ^l] . A com- 
bined analysis of both CMB and galaxy clustering is 
indeed necessary in order to break some of the degen- 
eracies between the cosmological parameters that are 
affecting the single datasets. Recent combined analy- 
sis, for example, have provided important constraints 
on the amplitude of the matter fluctuations (l5[ 
on the gravity waves background p0| ] , on the dark en- 
ergy equation of state ^3[ Q and on the neutrino 
mass [£5[ p6f . 

Thank to this new evidence, it is now accepted 
that the structure in the Universe is caused by an 
almost scale-independent primordial density pertur- 
bation, existing before cosmological scales start to en- 
ter the horizon. Nearly exponential inflation can eas- 
ily explain the origin of this perturbation, because it 
converts the quantum fluctuation of every light scalar 
field into a scale-independent classical perturbation. 
At the present time, there is no equally successful al- 
ternative explanation for the origin of this perturba- 



tion. In particular, there is no reason to expect the 
quantum-to-classical conversion to take place in extra 
space dimensions, though of course a string-theoretic 
description of such dimensions may play a vital in- 
direct role by determining the form of the effective 
four-dimensional field theory. 

In this paper, we combine the latest CMB and 
galaxy clustering data in order to place constraints 
on the inflationary model responsible for the primor- 
dial density perturbation. We will in fact assume 
the simple paradigm of a slowly-rolling inflaton field, 
whose perturbations are the seed for the primordial 
density perturbations |?7| |2£| 1 . If this hypothe- 
sis is correct, observation will be able to distinguish 
between the few simple and well-motivated inflation 
models that exist, because they give different predic- 
tions for the spectral index of the primordial den- 
sity perturbation. From among such models, we fo- 
cus in this paper on the running-mass model of in- 
flation 
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|32 , |33| , P4 |35| , which is the only one 
which can that can give a strongly scale-dependent 
(running) spectral index. Comparing with a suite of 
observations, we find that significant running is still 
allowed. 

A similar attempt in constraining this kind of scale- 
dependence was made in and 1 37 , where however 
only a limited subset of the CMB and LSS data then 
available was used. A model-independent approach is 
used instead in ]38| ] , where constraints in the plane n 
and its first derivative dn/dlak were obtained. The 
main difference with our approach is that there it is 
implicitly assumed that the scale dependence is small, 
while this is not necessary the case for the running 
mass model. 

Our paper is organized as follows: In section II we 
discuss the running mass model. In section III we 
present our data analysis method and results. Finally, 
in section IV, we discuss our conclusions. 



1 The alternative [Q] is to make a 'curvaton' field responsible, 
different from the inflaton if indeed the latter exists. 
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II. THE RUNNING MASS MODEL. 



A. The potential 

The running-mass model @ H || f§ || f§ 
seems to be the only one which is both very well- 
motivated from the particle physics point of view and 
also presents a strong scale-dependence of the spec- 
tral index. The model is based on supersymmetry 
and identifies the inflaton field with a flat direction 
of the supersymmetric scalar potential. Global su- 
persymmetry is broken during inflation by the large 
vacuum energy that drives inflation and so soft mass 
terms are generated for all scalar fields |59|. Taking 
the one-loop corrections to the tree-level potential into 
account, the most general expression for the potential 
along a flat direction, like the inflaton's, can then be 
cast into the form 



f 



V = V o + -m 2 (ln0)(f> 2 + 



(1) 



where m 2 is the soft inflaton mass and the dots in- 
dicate non-renormalizable terms which are supposed 
to be negligible because cf>/Mp is exponentially small 
(the usually-dominant renormalizable quartic term is 
absent precisely because we are using a flat direction 
in the scalar field space). The mass dependence on 
the renormalization scale, in our case identified with 
the value of the inflaton field, is given by the renor- 
malization group equation pC] 



dm 2 
d In c 



dm 2 
dlnQ 



(2) 



and therefore is proportional to the inflaton couplings, 
as given below in Eq. (||). 

Ignoring for the moment the running, a generic su- 
pergravity theory will generate a soft mass-squared 
with magnitude at least of order |m 2 | ~ Vo/Mp, which 
is marginally too big to support inflation. This is a 
problem for any model of inflation, whether or not the 
mass term is the one that is supposed to dominate. In 
all models except the running mass model, the prob- 
lem is solved either by imposing a global symmetry to 
protect the flatness of the potential, or by supposing 
that the mass-squared is accidentally suppressed. The 
running mass model instead accepts the generic value 
of the mass-squared at the Planck scale (or some other 
high scale), and relies on the loop correction to suf- 
ficiently reduce it in the regime where inflation takes 
place. 

For a particle with gauge and Yukawa interactions, 
we have that at one loop j3 m is given by Q 



2C ~ 2 D 
Pm = am + — — - A 
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m 



loop 



(3) 



where the first term arises from gauge particles loops 
and the second from matter loops. Above C, D are a 



positive group-theoretic numbers of order one, count- 
ing the degrees of freedom present in the loop, a is 
the gauge coupling, and fh is the gaugino mass, while 
A denotes a common Yukawa coupling and w 2 oop the 
common susy breaking mass-squared of the scalar par- 
ticles interacting with the inflaton via Yukawa inter- 
action. Note that the first term in Eq. (||) is always 
negative, while the second has no definite sign, since 
m ioop ' s defined as the mass squared splitting between 
scalar and fermionic superpartners and can have ei- 
ther sign. Also the case of a very-weakly-interacting 
inflaton gives (3 m — > and so that the constant mass 
potential is recovered. 

Over a sufficiently small range of </>, or for very small 
inflaton couplings, it is a good approximation to take 
a Taylor expansion of the running mass m 2 (In <fi) and 
then we have: 



V = V + ~m 2 (ln 



-'-^(M^M^), (4) 



where we have expanded around </>, and rescaled the 
last term w.r.t. Vo/Mp for future convenience. The 
dimensionless constant c is proportional to the mass 
beta function, 



M 2 f3 m /V 



(5) 



It has been shown |33| that for small c, as is re- 
quired by slow roll conditions, the linear approxima- 
tion is very good over the range of (j> corresponding to 
horizon exit for scales between /ccobe and 8/i -1 Mpc. 
In order to obtain also a crude estimation of the 
reionization epoch via the Press-Schechter formula, 
which involves the power spectrum at a scale of or- 
der k~ ei ~ 10~ 2 Mpc (enclosing the relevant mass of 
order 10 6 M Q ), we shall assume that the linear approx- 
imation is adequate down to this 'reionization scale'. 

For studying the inflaton potential, it is very useful 
to introduce a new parameter via 
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Then Eq. ^ takes the simple form p3[ 



V = V - --^W 2 I In — - - 

2M 2 v \ d>* 2 



leading to 
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(6) 



(7) 



(8) 



In typical cases the linear approximation is valid at 
4> — 4>*i and that point is then a maximum or a mini- 
mum of the potential. 

The running-mass model supposes that the relevant 
soft masses at the Planck scale (or some other high 
scale) have magnitude roughly of order 



I soft mass I ~ V 2 /M P - H , 



(9) 
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and that the relevant couplings (gauge or Yukawa) 
are very roughly of order 1. In general, one of the 
loops will dominate over the others making only one 
soft mass relevant besides the inflaton mass. With 
these assumptions, the coupling c defined by Eq. (S) 
is roughly 



lCT 1 to 10" 



(10) 



A bigger value of |c| is unviable because it would not 
allow inflation. On the other hand, using Eq. (^) as a 
very crude estimate of the inflaton mass at the Planck 
scale, one can see that a much smaller value of |c| is 
probably not viable either, since it would require the 
inflaton mass at that scale to be suppressed below the 
estimate Eq. (j|). 

For scalar masses, the estimate Eq. @ is expected 
to be valid provided that the following statements are 
true: (i) there is at most gravitational-strength cou- 
pling between the inflaton sector and the sector in 
which SUSY is spontaneously broken (ii) this break- 
ing comes from an F term as opposed to a D term 
(iii) during inflation (in contrast to the situation in 
the vacuum) the F term in the potential is not ac- 
curately canceled by the other, negative contribution 
to the potential. The third of these requirements is 
reasonable provided that the potential satisfies 



Vo > M 



s j 



(11) 



where Ms is the SUSY-breaking scale in the vacuum. 

The most economical assumption is that actually 
the mechanism of spontaneous SUSY breaking dur- 
ing inflation is basically the same as it is in the vac- 
uum, making Vo roughly of order Mg. For respectively 
the cases of anomaly-mediated, gravity-mediated and 
(low-energy) gauge-mediated transmission of SUSY- 
breaking to the Standard Model sector, this leads to 
the rough estimates 
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10 12 GeV 
10 10 GeV 
10 6 GeV . 



(12) 
(13) 
(14) 



The second case was assumed in the original works 
|3C| , |3S[| . It gives the correct normalization for 
the spectrum for reasonable parameters, but there is 
probably enough latitude to allow the first or even the 
third case. 

We emphasize that the estimate Eq. (JoJ) for the 
scalar soft masses in the inflaton sector is reason- 
able in all three cases, even though the soft masses in 
the Standard Model sector (evaluated of course in the 
vacuum, and of order 100 GeV) satisfy that estimate 
only in the second, gravity-mediated case. In the first 
(anomaly-mediated) case, this is because no-scale su- 
pergravity does not automatically suppress soft scalar 
masses during inflation, as it does in the vacuum ]27| . 
In the third (gauge-mediated) case it is because one 



can easily suppose that the gauge-mediation mecha- 
nism simply does not operate in the inflaton sector. If 
the relevant soft mass is a gaugino mass, the expec- 
tation of Eq. (0) is not so automatic because gaugino 
masses can be very small with some types of SUSY 
breaking, but it is not unreasonable either. 

In summary, we expect c in the range (10) in a 
typical running mass model of inflation. As we now 
see, this typically leads to significant scale-dependence 
of the spectral index. 



B. The spectrum and the spectral index 

Now we come to the predicted spectrum and spec- 
tral index of the primordial density perturbation. We 
suppose that it is generated by the inflaton field per- 
turbation, which means that it is purely adiabatic 
and gaussian. It is therefore specified by the curva- 
ture perturbation TZ(k), with k as usual the comov- 
ing wave number. This quantity is gaussian and hence 
specified by its spectrum Vn{k). 

To give the prediction of the running mass model 
for the spectrum, it is convenient to define 



cln(</>*/</>, 



COBEj 



(15) 



where the subscript COBE denotes the epoch of hori- 
zon exit for the COBE scale fc C OBE ^ 0.002(/i/ Mpc). 
Instead of the spectrum Vnik) we work with the 

At the COBE 



conventional quantity Sh = f^)!/ 2 - 
scale, the prediction of the running-mass model is 



Sn(kcOBE) 



1 



-\a\- l exp(a/c). (16) 



As already mentioned, the observed value 
Sh (^cobe) = 1.9 x 10~ 5 is easily obtained, for 
choices of Vo and that correspond to reasonable 
particle physics assumptions. 

In this paper we concern ourselves with the dra- 
matic scale-dependence of the spectrum, which is 
given by 
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2s / 

TV 



cAN(k) 



1 



2cAN(k) 



(17) 

where AN(k) = V CO be - N(k) = ln(fc/fc C OB E )- We 
plot the primordial power spectrum for the running 
mass model in comparison to the constant spectral 
index case in Figure [l| 

To discuss the spectral index and its running, we 
need the first four derivatives of the potential, 
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(18) 
(19) 
(20) 
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FIG. 1: Primordial power spectrum on a logarithmic plot 
for the case of a scale invariant spectral index n = 0.95 
(solid line) and n = 1.1 (dashed-dotted line), and for the 
running mass prediction with iicobe = 0.95 and two dif- 
ferent Wcobei i- e - n coBE — 0.04, corresponding to c = 
0.1545, s = 0.1295 (upper dashed line) or c = -0.1295, s = 
—0.1545 (lower dashed line), and ticobe = 0.01, corre- 
sponding to c = 0.08431, s = 0.05931 (upper dotted line) 
or c = -0.05931, s = -0.08431 (lower dotted line). 
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(21) 



Introducing s, the first four flatness parameters 
p8[ are (setting V = Vq in the denominators) 



V 



,cAN(k) 
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(22) 

(23) 
(24) 
(25) 



The parameters are evaluated at the epoch of horizon 
exit for the scale k. The first parameter e is negligible 
because (fr/Mp is taken to be exponentially small. The 
condition for slow-roll inflation is therefore just <C 
1, which is satisfied in the regime <fi ~ 4>* provided that 
|c| <C 1. This corresponds also to requiring |s| <C 1. 

Additional and generally stronger constraints on s 
follow from the reasonable assumptions that the mass 
continues to run to the end of slow-roll inflation, and 
that the linear approximation remains roughly valid. 
Discounting the possibility that the end of inflation 
is very fine-tuned, to occur close to the maximum or 
minimum of the potential, we have the lower bound 



s >e 



-cNcoe 



(26) 



Note that for negative c, this constraint is very strong, 
requiring a very large value of s even for small c and a 



FIG. 2: Theoretical expected valued of the parameters c, s 
for iVcoBE = 50; the solid-line-hatched region is strongly 
excluded by naturality assumptions, while the dashed-line- 
hatched region is only weakly excluded. The dotted line 
shows the prediction for the simplest case, s — e _cJVcOBE , 
where the linear approximation is valid up to the end of 
inflation, triggered by r\ — 1. 



kind of fine-tuning between s and c to give a reason- 
able value of n — 1 . 

For positive c, we also obtain a significant upper 
bound by setting AN = iVcoBE in Eq. (|30|), and re- 
membering that slow-roll requires |n — 1| < 1: 



s <e 



-cJVcOE 



(c> 0) . 



(27) 



In the simplest case, that slow-roll inflation ends when 
ri — l actually becomes of order 1 , this bound becomes 
an actual estimate, |s| ~ 6 -cN C obe_ j± s discussed in 
p(| , this upper bound can be relaxed for positive s if 
the running of the mass ceases before the end of slow- 
roll inflation. The approximate region of the s versus 
c plane excluded by these considerations is shown in 
Figure |2|. 

Since e is negligible, the spectral index is 



This gives J3C 



i(k) = 



i(k) - 1 



1 + 2?y - 6e 
I + 277. 



se cAN{k)_ c _ 



(28) 
(29) 



(30) 



The first derivative of the spectral index is given by 

dn(fc) =2sce cAiv( fc )_ (31) 
dhiK 

Clearly the spectral index is not constant within cos- 
mological unless s or c is very close to zero. 
At the COBE scale, 



n-coBE - 1 = 2(s - c) 
= 2sc. 



'COBE 



(32) 
(33) 
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The straight line s = c in the s vs c plane corresponds 
to ncoBE = 1< The scale independent case of constant 
n = 1 is given by the origin s = c = 0, while constant 
spectral index different from 1 is realized either on the 
c = axis for s — (n— l)/2 or on the s = axis for 
c= -(n- l)/2. 

Taking cosmological scales to span a range AN ~ 
10, we see that the running mass model generates a 
change Sn ~ lOsc in the spectral index. This change 
should eventually be detectable unless both \c\ and \s\ 
are towards the bottom of their expected ranges. 

We will show in the following the allowed region 
both in the s, c parameter space and in the ^cobe vs 
"■cobe — 1 plane. Note in this respect that not all val- 
ues of ^cobe are allowed in the running mass models: 
from the definition above, we obtain the constraint 



l COBE 



> - 



(riCOBE _ !) 2 



(34) 



so that a decreasing spectral index ispossible only if 
"-COBE is different from 1. Also Eq. ( |33|) is symmetric 
under reflection with respect to the s + c = line, so 



fixing the variables ticobe — 1 and 



^COBE 



determines 



two values of s,c. We will investigate in the following 
how strongly the data are sensitive to a variation of 
n' and therefore able to disentangle this degeneracy. 
From all this, it follows that fitting for arbitrary value 
of n' is not exactly equivalent to performing a fit for 
the running mass models. 

These predictions for the spectral index and its 
derivative are to leading order in slow-roll. To second 
order, the spectral index is given to a good approxi- 
mation by p7[ 



n(k) - 1 



s(l + 1.06c)e 



cAN(k) 



(35) 



We neglected a term of order rj 2 , and the remaining 
(unknown) error is expected to be of order tr 3 . Both 
of these are negligible which means that the second- 
order correction should be rather accurate. We shall 
not include it though, because it just corresponds to a 
small rescaling of the quantity s, whose precise value 
depends on unknown physics p7[ . Note that it is in 
any case easy to obtain the results at next-to-leading 
order, by rescaling s by the factor 1/(1 + 1.06c) or 
shift n C OBE - 1 by 1.06n' COBE . 



III. OBSERVATIONAL CONSTRAINTS 



A. Method 

We compare the recent cosmological observations 
with a grid of theoretical models computed with 
the public available CMBFAST |Sj code. We re- 
strict our analysis to flat, adiabatic, A-CDM mod- 
els and we sample the various parameters as follows: 



fl cdm h 2 = uj cdm = 0.01, ...0.40, in steps of 0.01; 
tt b h 2 ee uj b = 0.001, ...,0.040, in steps of 0.001 and 
tt A = 0.0, ...,0.95, in steps of 0.05 The va lue of the 

Hubble constant is given by h = "^"'^^ 

We allow for a reionization of the intergalactic 
medium by varying also the compton optical depth 
parameter r c in the range r c = 0.0, 0.20 in steps of 
0.05. Greater values of r c are in disagreement with re- 
cent estimates of the redshift of reionization zr ~ 6±1 
(see e.g. ^J) which points towards r c ~ 0.05. As dis- 
cussed later, given a cosmological model, we checked 
for the consistency of the assumed optical depth with 
its expected values derived from a simple reionization 
scenario. 

We vary the 2 parameters of the running mass 
model in the range -0.30 < c < 0.30 and -0.47 < 
s < 0.47 in a 50 x 50 grid. To save computing time, 
we compute the CMB anisotropies transfer functions 
for each cosmological model just one time and then 
we integrate them different times looping over c and 
s. 

For the CMB data, we use the recent results from 
the BOOMERanG-98, DASI, MAXIMA-1,CBI, and 
VSA experiments. The power spectra from these ex- 
periments were estimated in 19, 9, 13, 14 and 10 bins 
respectively (for the CBI, we use the data from the 
MOSAIC configuration), spanning the range 2 < I < 
3500. However, since in this work we are interested 
only in comparing the data with the expected pri- 
mary anisotropies, we limit our analysis to the region 
I < 1500 which is likely not to be affected by sec- 
ondary effects like Sunyaev-Zel'dovich (see e.g. Q). 

The likelihood for a given theoretical model is de- 
fined by -21n£ = {Cf-Cf )M B B'{Cf, -Cf?) where 
Mbb' is the Gaussian curvature of the likelihood ma- 
trix at the peak and Cb are the experimental (theo- 
retical) band powers. We discard the first bin of the 
CBI dataset (0 < £ < 400), due to the asymmetric 
window function and the high sample variance. 

We consider 10%, 4%, 5%, 3.5% and 5% Gaussian 
distributed calibration errors for the BOOMERanG- 
98, DASI, MAXIMA- 1, VSA, and CBI experiments 
respectively and we include the beam uncertainties 
by the analytical marginalization method presented 
in @. 

In addition to the CMB data we incorporate the 
real-space power spectrum of galaxies in the 2dF 100k 
galaxy redshift survey using the data and window 
functions of the analysis of Tcgmark et al. [Q . 

To compute C 2dF , we evaluate pi — P(ki), where 
P(k) is the theoretical matter power spectrum and 
ki are the 49 k- values of the measurements in Ejj. 
Therefore we have -2lnC 2dF = £JP;- (Wp)i\ 2 /dFf, 
where Pj and dPi are the measurements and corre- 
sponding error bars and W is the reported 27 x 49 
window matrix. We restrict the analysis to a range 
of scales where the fluctuations are assumed to be in 
the linear regime [k < 0.2/i _1 Mpc). When combin- 
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ing with the CMB data, we marginalize over a bias b 
considered to be an additional free parameter. 

We also include the recent 13 data points on the 
matter power spectrum obtained by Croft et al. ||45f 
using Lyman a forest data from 53 quasar spectra. 
The theoretical transfer functions are in this case com- 
puted up to redshift z = 2.72 and compared by a sim- 
ple chi-square analysis with the data, including the 
±25% error in the data amplitude. 

Moreover, following a similar approach used in a 
previous analysis |37j we use reionization bounds. 
Practically, we deduce for each theoretical model, the 
expected value of rf 1 , 



0.035 n b h 



(V(l-«A)(l + *Ji) 3 + nA-l) 



from the reionization redshift zr estimated using a 
Prcss-Schechter formula. Taking / to be the frac- 
tion of matter collapsed into objects with mass M = 
10 6 Af o , 1 + z R is given by 



zr 



V2a{M) 
fl (l-f2 A ) 



erfc- x (/) 



/ — 1 



(37) 



where a(M) is the present, linearly evolved, rms den- 
sity contrast with top- hat smoothing and S c = 1.7 
is the overdensity contrast required for gravitational 
collapse; g(\ — is the suppression factor of a(M) 
at the present epoch for the case of a non vanishing 
cosmological constant: 

g ai) = ln(- + —n - —n 2 + n 4 ^ 7 ) ' . (38) 

yv ' 2 \70 140 140 / V ' 

In the case of the running mass model, a(M) can 
become very large even for moderate values of c, 
if chi(kfi/kcoBE) — 1, so that the reionization con- 
straint is very powerful. To give a feeling of this strong 
dependence, we show in Figure ^| the value of Zr as a 
function of c for the simple cases s — c, s = c — 0.05 
at fixed cosmological parameters. 

For each model in the database with a given op- 
tical depth r c we then consider the further prior 
2lnC rel = [t c -t c 4/i ] 2 /0.05 2 where the 1-cr error should 
take into account the uncertainties due to the different 
reionization mechanisms. 

Finally, we also include a constraint on the variance 
a on a mass scale of 1O 6 /i _1 M of 9 ± 1 consistent 
with the limits on the linear power spectrum obtained 
using measurements of substructure in gravitational 
lens galaxies in the analysis of |4(| . 

One should keep in mind that there are many 
caveats in the inclusion of all above constraints. 

We attribute a likelihood to each value of c and s by 
marginalizing over the nuisance parameters. We then 
define our 68% (95%), confidence levels to be where 
the integral of the likelihood is 0.16 (0.025) and 0.84 
(0.975) of the total value. 
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FIG. 3: Theoretical expected valued of zr for / ~ 1 as a 
function of c for the case s = c (dashed curve) and s = 
c — 0.05 (dot-dashed curve). The cosmological parameters 
are chosen as h — 0.72, ui c dm = 0.086 and cot = 0.02. We 
show also the reference line zr = 6 (solid line). 



B. Results 



The likelihood contours in the c — s plane obtained 
analyzing the CMB data under the assumption of a 
set of "weak" priors (h — 0.65 ±0.2, to > HGyrs) are 
plotted in Figure ^. As we can see, even if the space 
of models analyzed is quite broad, the CMB data is 
able to give strong constraints along the s — c direc- 
tion. In particular, we obtain the constraint ncoBE = 
1 + 2(s - c) = 0.961^ and n$ = at 68% c.l... We 
found that including stronger priors {h = 0.72 ± 0.08, 
Sl m — 0.3 ± 0.1) on the nuisance parameters does not 
change these results and does not improve significantly 
the constraints on ticobe- Assuming negligible reion- 
ization (r = 0) gives ucobe = 0.93lo;o3- 

Also plotted in Figure || (Top Panel) are the like- 
lihood contours obtained in the combined 2dFGRS- 
CMB analysis. As we can see, even the inclusion of 
the 2dF data does not improve significantly the CMB 
constraints. The combination of the 2dF data is ex- 
tremely powerful in constraining parameters like f2 m 
and h, that define the epoch of equality and the po- 
sition in the k space of the matter spectrum. How- 
ever, in the CMB+2dF combined analysis, we found 
ucobe = 0.97i°;EJl at 68% c.l.. 

In the center panel of Figure [| we report the likeli- 
hood contours obtained in the combined CMB+Ly-a 
analysis. As already noticed by [[38], we find that the 
Lyman a data are able to restrict more strongly the 
scale dependence of the spectral index and therefore 
exclude the parameter space at large |c|; in particular 
at 68% we have sc < 0.026, so that a variation of the 
spectral index over cosmological scale of the order of 
10% is allowed. 

Finally, in the bottom panel of Figure [| we plot 
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FIG. 4: Likelihood contours in the c — s plane. The 
gray contours are the 95% c.l. from the CMB anal- 
ysis. The dark contours are from CMB+2dF analy- 
sis (Top), CMB+Lya (Center Panel), CMB+Reionization 
constraint (Bottom Panel). The straight lines corresponds 
to ncoBE = 0.8,1.0 and 1.2. 
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FIG. 5: Likelihood contours in the (ucobe — 1) — (^cobe) 
plane for different datasets. 



the likelihood contours obtained in the combined 
CMB+reionization analysis. Also in this case part 
of the region of large s,c is excluded, but still a siz- 
able variation and values of |c| up to 0.2 are allowed. 
It should be emphasized that this reionization con- 
straint is only one of theoretical consistency. As one 
can see from Figure [|, the self-consistent value of Zr 
becomes very large for c > 0.10. If it turns out that 
Z R % 10) the top part of the allowed region 'CMB + 
Reionization' will be excluded. 

Focussing on the region (s + c > , c > 0), we 
plot the likelihood contours in the 2(s — c) vs. 2sc 
plane in Figure ^. As we explained before, 2(s — c) 
gives the value of the spectral index ucobe — 1 on 
COBE scales, while 2sc = n' COBE gives the bend in 
the spectrum. 

As we can see from the figure, the recent CMB data 
alone constrains «cobe < 0-1 a ^ 95% C.L., while in- 
cluding the constraints from 2dF, Ly-a, reionization 
redshift, limits the amount of deviation from scale in- 
variance to n' COBE < 0.05 at 95% c.l.. Adding the 
strong lensing constraint from [^6[ , as discussed ear- 
lier, further constrains n 1 < 0.04 at 95% c.l.. It 
is important to notice that a lower spectral index 
n < 0.95 is consistent with the lensing observation 
only if «cobe > 0- Note also that in all cases the best 
fit value of n' COBE is not in the centre of the allowed 
region, but lower, towards the value ^cobe = 

0. Fur- 
thermore, the inclusion of a bend in the power spec- 
trum, does not seem to affect in a considerable way the 
CMB constraint on n since there is a weak correlation 
between these two variables. 

The correlation between n' C Q BE and the remaining 
parameters is investigated in Figure ^ where the like- 
lihood function for n' COBE is plotted in function of cob, 
oJcdm and the Hubble parameter, h. 

As we can see, the correlation with these parame- 
ters is extremely small, meaning that the inclusion of 
a bending in the primordial spectrum would not dras- 
tically affect the actual constraints on the parameters 
considered. However, an increase in ^cobe can a ^ 
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FIG. 6: Correlations between riQ OBE with different pa- 
rameters. From top to bottom, we consider the physical 
density in baryons u>b, in cold dark matter ui c dm and the 
Hubble parameter, h. 



low larger (smaller) values of cjf, (us c d m )- The Hubble 
parameter h is poorly constrained by CMB alone. 

We have also considered the effect of variations in 
ncoBE and n coBE 011 the predicted level of rms mass 
fluctuations in spheres of 8/i _1 Mpc. At the moment, 
there is some sort of tension between the possible val- 
ues of this parameter obtained from local cluster abun- 
dances (erg ~ 0.75 for Qm — 0.3, see e.g. Wh) and 
weak lensing (cr% ~ 0.9 for Qm = 0.3 see e.g. |H8j|). 

In Figure we plot the likelihood contours from the 
CMB analysis for erg against ticobe — 1 (Top Panel) 
and ri.Q OBE (Bottom Panel) respectively. As expected, 
increasing both ticobe — 1 and "cobe increases the 
level of erg , but the correlation is not relevant enough 
to actually detect a bending in the spectrum with 




FIG. 7: Likelihood contours in the (ticobe — 1) — era (Top) 
and n' COBE — as planes. 



present observations. 

We are in qualitative agreement with other recent 
constraints on a running spectral index ]38) . Precise 
comparison is difficult mainly because we did not con- 
sider a background of gravity waves (that can weaken 
our final result) since in the running-mass models this 
is negligible. Furthermore, we considered a more up- 
dated CMB dataset respect to the one considered in 
|3S|| and we considered different datasets other than 
CMB. Finally, j3£| assumes that n' can take any 
value, whereas the running-mass model excludes sig- 
nificantly negative n' (Eq. (|34|)). The running-mass 
model also allows variation of n'(k), but this is not 
very significant at the present level of observation. 

We end this discussion of our results by comment- 
ing on a recent analysis by Caprini, Hansen and Kunz 
p9|] which investigates the variation of n'(k) that is 
allowed by observation. First, they assume that the 
flatness parameter a 3 is constant to the end of slow- 
roll inflation, and taking A^obe = 50, they find 
ct 3 > -3.5 x 10~ 4 . The running-mass prediction 
Eq. (|25|) for cr 3 is roughly constant to the end of slow- 
roll inflation if |c| < 1/50, and in this regime one 
needs \s\ < 0.1 to have n sufficiently close to 1; com- 
bining these results gives roughly |er 3 | < 4 x 10 -4 
in agreement with pSfl . However, assume next 
that instead V"" is roughly constant, and derives a 
bound Mp/"" jV > -0.02. In contrast, the running- 
mass model predicts MpV""/V = c(M P /<j)) 2 which is 
exponentially large (and of either sign) and strongly 
varying. Thus the bound on V"" of Q does not ap- 
ply to the running-mass model in any regime of its 
parameter space. We remind the reader again that, 
at the moment, the running-mass model is the only 
well-motivated model giving significant running of the 
spectral index. 



IV. CONCLUSIONS 

The rather full analysis that we have described con- 
firms the general picture indicated by previous anal- 
yses [[37|. The allowed region in the c vs s plane de- 
picted in Figure ^ should be compared with the re- 
gion shown in Figure H which approximately delin- 
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eates the theoretically disfavoured region, and also 
with the minimum value |c| > 10~ 2 which is probably 
needed to generate enough running of the mass even 
if we go from the Planck scale to 100 GeV. Combining 
all of these, we see that if |c| is significantly above the 
minimum value, only the version of the model with c 
and s both positive is viable. In that case, the spectral 
index has significant running which will be detectable 
in the forseeable future. On the other hand, if |c| is 
really of order 10~ 2 , all choices of the signs of c and 
s are possible except maybe negative c with positive 
s. Furthermore, if that extreme case can be realized 
in a viable running-mass model the running of n will 
be so small that it may never be detectable. 

Looking at the observational situation in more de- 
tail, our results show that the CMB data can put very 
strong constraints on the value of the spectral index 
at large scales, ucobe = l + 2(s — c), but still allow a 
pretty large scale-dependence. Other information on 



the power spectrum, like Lyman a data, or strongest 
assumptions about the reionization epoch are needed 
to reduce the parameter space in the s + c direction. 
Still values of |c| of the order of 0.2, and therefore 
n coBE °f the order 0.04, are allowed. Note also that 
our allowed region is more or less symmetric under re- 
flection with respect to the s + c = line: this means 
that the present data are not sensitive enough to dis- 
tinguish the variation of n' that is predicted by the 
running-mass model. Finally, we show that the inclu- 
sion of the running in the spectral index has a mod- 
erate effect on the present CMB determination of the 
physical matter densities in baryons uib and cold dark 
matter Lo cdm - 
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